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We study the energy levels of non-interacting electrons confined to move in two-dimensional 
billiard regions and having a spin-dependent dynamics due to a finite Rashba spin splitting. The 
Green's function for such Rashba billiards is constructed analytically and used to find the area 
and perimeter contributions to the density of states, as well as the smooth counting function. We 
show that, in contrast to systems with spin-rotational invariance, Rashba billiards always possess a 
negative energy spectrum. A semi-classical analysis is presented to interpret the singular behavior 
of the density of states at certain negative energies. Our detailed analysis of the spin structure of 
Rashba billiards reveals a finite out-of-plane spin projection for electron eigenstates. 
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I. INTRODUCTION 

Spin-dependent phenomena in semiconductor nanos- 
tructures have attracted great current interest In- 
triguing effects can arise in non-magnetic systems due to 
the presence of spin-orbit coupling. Structural inversion 
asymmetry in semiconductor heterostructures has been 
shown^ to give rise to a spin splitting of the same type 
as was discussed in an early paper by Rashba^ Its tun- 
ability by external gate voltages^i has motivated the 
theoretical design of a spin-controlled field-effect tran- 
sistor. 8 Novel spin properties arise from the interplay 
between Rashba spin splitting and further confinement 
of two-dimensional electrons in quantum wires^iiSiiiiiS 
ringS) 13 i 14 or dotsi 15 i 16 ' 17 i 18 i 19 i 20 ' 21 Spin-orbit coupling 
has also been shown to affect the statistics of energy levels 
and eigenfunctions as well as current distributions^^ 
The interplay between spin-orbit coupling and external 
magnetic fields was analyzed theoretically using random 
matrix theory^* 

In this work, we study Rashba billiards, i.e., non- 
interacting ballistic electrons moving in finite two- 
dimensional (2D) regions whose dynamics is affected 
by the Rashba spin-orbit coupling. In the one-band 
effective-mass approximation, the Hamiltonian with 
Rashba splitting in 2D is given by2£ 



H = 

H = 
U = 



n 



P 2 y 



2m* 

<?xPy ~ VyPx 



(la) 

(lb) 
(1c) 



where a x ,a y are Pauli matrices. This Hamiltonian gov- 
erns the electron dynamics inside the billiard with Dirich- 
let boundary conditions at the perimeter (see Ref ll8|) . 
The Rashba spin-orbit coupling strength a can be con- 
veniently measured in terms of a wave-number scale 



k so = m*a/h 2 , The spin-precession length defined as 
Tr/k so can be tuned independently of the system size££*£ 
Furthermore, the tunability of the Rashba spin-orbit 
coupling strength is a convenient tool to induce changes 
of the billiard's energy spectrum without applying exter- 
nal magnetic fields. 

One of our central quantity of interest is the density 
of states (DOS) g(E) = J2 n ~ E n ) and the counting 
function N(E) = Q(E-E n ) for Rashba billiards with 
energy levels E n . Here 8{x) and 0(x) are the Dirac delta 
function and the Heaviside function, respectively. The 
density of states and the counting function of normal 
billiards (without spin-orbit coupling, i.e., for a = 0) 
have been extensively studied in the literature. They can 
be derived from the Green's function of the system. The 
smooth counting function N(E) is given by the so-called 
Weyl formulaj2£i2L2£ which is an asymptotic series of the 
exact counting function N(E), in terms of the energy E. 
A good introduction to this problem is given by Baltcs 
and Hilf, 2 £ and recently by Brack and Bhaduri^ and 
many applications can be found in Refs 29 3(ill3lll32ll33t 

a 

We have derived the first two leading terms in the Weyl 
formula of the smooth counting function N(E) for arbi- 
trary shapes of Rashba billiards. Our approach is based 
on the image method of Berry and MondragonS devel- 
oped for neutrino billiards, which have two-component 
wave functions and in this respect are rather similar to 
the Rashba billiards discussed here. We will show that 
the first term of N(E) is proportional to the area of the 
billiard, while the second one is proportional to the length 
of the perimeter of the billiard. Moreover, we find that 
the density of states is singular at the bottom of the 
spectrum. This singular behavior occurs independently 
of the billiard's shape and is most striking if the Rashba 
parameter is large. 

The circular Rashba billiards is the simplest of confined 



2 



systems that can be treated analytically pj&SLiS Follow- 
ing the approach outlined above, we also calculate the 
smooth counting function N(E) for circular Rashba bil- 
liards, and besides its first two leading terms (which co- 
incide with the results derived for arbitrary shapes of 
Rashba billiards) we give higher-order correction terms. 
Our analytical result for N(E) and that obtained from 
the numerically calculated exact energy levels are in per- 
fect agreement. 

In the absence of any lateral confinement, the energy 
dispersion for the Rashba Hamiltonian (JTJ splits into two 
branches^: 



E(k x , k y ) 



2m* 



(|k| ± fc so ) - k 2 s , 



(2) 



where k = (k x , k y ). The spin splitting is a consequence 
of broken spin-rotational invariance. The spin of energy 
eigenstates, which are labeled by a 2D vector k, is po- 
larized perpendicularly to k,— Hence, no common spin 
quantization axis for single-electron states can be defined 
in the presence of spin-orbit coupling. As can be seen, in 
the range < k < 2fc so , one branch has negative energies 
bounded from below by — A so = — h 2 k^ /(2m*). 

Similarly, a laterally confined 2D system in the pres- 
ence of Rashba type spin-orbit interactions has also a 
negative energy spectrum. In this paper, we present 
interesting features of the energy spectrum for circular 
Rashba billiards, focusing especially on its negative en- 
ergy eigenvalues. We have found that for a circular shape, 
the density of states has additional singularities at neg- 
ative energies. We obtain analytic results for their posi- 
tions. Their corresponding eigenspinors have a finite spin 
projection in the direction perpendicular to the billiard 
plane, which is the direct result of imposing hard- wall 
boundary conditions. 

Results presented in this article extend work reported 
in Ref. 39. Its organisation is as follows. The proper- 
ties of arbitrarily shaped Rashba billiards are discussed 
in Sec. [H] We present an algebraic expression for the 
free-space Green's function in the presence of Rashba 
spin-orbit coupling in Sec. Ill Al Susbequently, in subsec- 
tion ^TO] the first two leading terms of the Weyl formula 
are derived. Using the eigenstates presented in Sec. Ill CI 
in the absence of lateral confinement, we derive an alter- 
native expression for the free-space Green's function in 
polar coordinates in Sec. Ill Dl The circular Rashba bil- 
liards are discussed in Sec. lIIII An analytical formula for 
the Green's function is derived in Sec. lIII Al for this case, 
while the derivation of the smooth counting function is 
presented in Sec. lIII 51 including its comparison with the 
numerically calculated result. For negative energies the 
counting function is calculated in Sec. IIII CI while the 
spin structures is discussed in Sec. IIII Dl Finally, our 
results are summarized and conclusions given in Sec. II VI 



II. ARBITRARY SHAPES OF RASHBA 
BILLIARDS 

In this section we derive the smooth part of the density 
of states and the smooth part of the counting function, 
i.e., the two leading terms in the Weyl formula 26 - 27,28 
for arbitrary shapes of Rashba billiards. These smooth 
functions are obtained by averaging the exact DOS and 
counting function over a small energy range around an 
energy E. 

The exact density of states g(E) expressed in terms of 
the retarded Green's function (see e.g.,)2& is given by 

q(E) = — lim ImTrG(£ + ij),r,r'), (3) 

where the trace means the limit r — * r', integration of 
r over the area of the billiard, and the trace in spin 
space. The exact Green's function G(z, r, r') is the 
position representation of the Green operator G(z) = 

(z — H) , which in addition, satisfies the boundary con- 
ditions. Then, the exact counting function is defined by 
N(E) = J^ oo g(E')dE'. 

Usually, the exact Green's function satisfying the 
boundary conditions is not known. However, one can al- 
ways write the exact Green's function as a sum of the so- 
called free-space Green's function and a correction with 
which the exact Green's function satisfies the boundary 
conditions. The free-space Green's function G<x,(E, r, r') 
is the Green's function of the infinite system and does 
not satisfy the boundary conditions at the boundary of 
the billiards. In this paper, we calculate the free-space 
Green's function G oa (E, r, r') for the Rashba Hamilto- 
nian and in case of circular Rashba billiards, the ex- 
act Green's function which satisfies the Dirichlet bound- 
ary conditions at the boundary of the billiards. 

The first term in the Weyl formula, called area term, 
can be obtained by replacing the exact Green's func- 
tion with the free-space Green's function G OQ (£', r, r') in 
Eq. It is always proportional to the area of the bil- 
liard. Higher-order terms in the Weyl formula are the 
corrections to the area term taking into account the exact 
Green's function. The smooth part of the first correction 
term is called perimeter term because it is proportional 
to the length of the perimeter of the billiard. 



A. Free-space Green's function for Rashba billiards 

All our subsequent calculations are crucially based 
on the knowledge of the free-space Green's function 
G aD (E, r, r') for Rashba billiards. In this subsection, we 
we present its derivation. 

At a given energy E, two propagating modes exist 
whose wave vectors can be found from the dispersion re- 
lation (01 : 



|k| = k± = \k T kso\, where k 



2m*E 

^^ + k l- (4) 
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Using the identities for the Pauli matrices one can eas- 
y 



ily show that U 2 = p 2 + p 2 and the Rashba Hamiltonian 
can be written as 



2m* h 



The free-space Green operator (E) = (E — H) cor- 
responding to the Rashba Hamiltonian reads then 



The retarded Green's function in position representa- 
tion is given by 



G 00 (E,r,r') = {r\G 00 (E + ir ] )\r'), (9) 



where E is a real number and 77 — > + . The two terms in 
Eq. JSJ in position representation involve two functions: 



Goo{E) 



2m* 



k 2 (E)-l^ + k 



(5) 



where k(E) is given by Eq. Here E can be a complex 
number. Using the operator identity 



A - A' 



2\ 



x - n 1 + ( a - i N ~ ' 



. (6) 



where A is a scalar and A is an operator, one can decom- 
pose G 00(E) as 



dco[E) = w 



[sgn(E)k, 



, (7) 



where k± are given by Eq. (JIJ. Now using the operator 
identity (a± A\ * = (a =p A\ (\ 2 - A 2 ) , one finds 



Goo (E) - t , 



m* 1 



U 



k 2 - 



Sg a(E)k + + ^\ (k\-t^ 



(8) 



kl-Zp) \r'). (10) 



After a simple limiting procedure one can show that 

k 2 + {E + iri) = k 2 + (E)+sgn{E)im (11a) 
k 2 _{E + iri) = k 2 _(E)+ir). (lib) 

The two functions in (|lf)|> can be evaluated by the fol- 
lowing identities (see e.g.,) 2 ^: 



r| ( k 2 - ^ ± it] J |r' 



-i^(fe|r-r'|), 
-^ 2) (fc|r-r'|), 



(12) 



(1 2) 

where Hq ' (x) are the Hankel functions of zero order, 
and k > 0. 

Finally, using Eqs. J8J - (|12f> we can easily find 



n it? <\ - im * 1 
G x (E,r,r) = - 1T - 



'(k- - |) >(*-|r- r'|) + (fc+ + |) H^(k + \r- r'|)] , for £ > 0, 
" - f ) ^(fc-lr - r'|) - (-*+ + I ) ^ 2) (fc + |r - r'|)l , for £? < 0. 



(13) 



We note that, for negative energies E, the retarded 
Green's function contains incoming circular waves be- 
sides outgoing waves. 



tor U is zero since U is an off-diagonal matrix in the spin 
space. Then, it is easy to see that the leading term in 
the DOS becomes 



B. Area and perimeter terms of the density of 2m* 1 \ k, for E > 0, 

States 6>arca(^) = " To" 7 { (14) 

2n R k { k so , for £ < 0, 

First, consider the area term of the Weyl for- 
mula. Now, in Eq. @ we replace the exact Green's 

function G(E, r, r') by the free-space Green's function where A is the area of the Rashba billiard. Therefore, 

Goo(E, r, r') given by Eq. (|13|) . The trace of the opera- the integration of the DOS yields the counting function: 
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_ A 2m* 

7T h 




for E > 0, 



for - A so <£:< 0. 

(15) 

It follows directly from Eq. (|15|l that, for negative ener- 
gies, the DOS shows a 1/ y/E + A so singularity at the bot- 
tom of the spectrum E — > — A so . The area term (I15|l can 
alternatively be derived from the classical phase-space 
integral in the underlying classical approach. However, 
the classical dynamics of electrons in Rashba billiards is 
described by two Hamiltonians^ which are reminiscent 
of the two dispersion branches The constant-energy 
surfaces in phase space are different for the two Hamil- 
tonians, yielding different contributions to the classical 
phase-space integral. This simple calculation also leads 
to Eq. 

For arbitrary shapes of Rashba billiards, we can also 
determine the perimeter term of the DOS and the count- 
ing function. This term can be derived from the general- 
ization of the image method of Refl30l using only the free 
space Green's function. The actual calculation is very 
much similar to that applied by Berry and Mondragon 35 
for neutrino billiards. The Dirichlet boundary conditions 
can be approximately satisfied by regarding the boundary 
as straight and using the approximate Green's function 

G(r,r') « G ot) (r,r') + G<(r,r') 1 where (16a) 
G,(r,r') = -Goofari), (16b) 

and Ti is the image of r' on the boundary outside the 
billiard. Obviously G(r, r') is still a solution of the 
Schrodinger equation in the variable r. To calculate 
the trace in J3J of Gj(r,r') we define r = (n, s) and 
r' = (n, s+cr), where n and s are the coordinates of r per- 
pendicular to and along the boundary. Of course n < 



trace corresponds to a — ► 0. Now, | r— Vi |= y (2n) + a 2 
and the correction to the DOS, i.e., gperim(E) coming 
from Gi(r, r') can be written as 

Qpci-im(E) = -- lim ImTr G t (E + ir/,r,r') 

7T 7)^0+ 



1 TO* 1 



ds 



dn 



2tt h 2 k Jo 
[k-J (k-2n) + k + J (k + 2n)}, (17) 

where the factor 2 of the trace in the spin space has 
already been included. Using the integral J °° Jo(ax)dx = 
1/a with a > we obtain 



£*perim(-^) 



L 2 



in 



1 

^k' 



(18) 



4tt hf- 

valid for all energies E > — A so . Here £ is the length of 
the perimeter of the billiard. Finally, the integration of 
the DOS yields the counting function: 



N n 



(19) 



valid also for all energies E > — A so . The minus sign is 
a consequence of Dirichlet boundary conditions. 

In summary, the first two terms in the Weyl for- 
mula for arbitrary shapes of Rashba billiards reads as 
N(E) = N aTca (E) + N pelim (E). Note that for zero spin- 
orbit coupling, N(E) coincides with the previously de- 
rived result for 2D billiardo 27 i 28 i 29 i 31 i 32 (apart from a fac- 
tor 2 due to spin). 



C. Eigenstates for infinite systems 

In this section the eigenvalues and eigenstates of the 
free-particle Rashba Hamiltonian given by Eq. JIJ are cal- 
culated in polar coordinates. These results will be used 
in the subsection 111 Dl to rewrite the free-space Green's 
function [Eft in a form which is suitable for calculations 
in case of circular Rashba billiards presented in Sec. IIIII 

The Hamiltonian JQ) can be rewritten in polar coor- 
dinates r = (r, ip) and we have H = Hq + ? U, where 



U 
h 



Id Id 2 



2m* \ dr 2 r dr r 2 dtp 2 



- ' ' #■ 

or 



i_d_ 

r dip 



-up ( _d i_ d \ 

V dr r dip J 



(20a) 



I. (20b) 



Since the Hamiltonian H commutes with the total angu- 



lar momentum operator J 2 



-ihdu: 



I a z , the station- 



ary Schrodinger equation H \ x) = E \ x) can be solved 
using the following ansatai&iS 



GiZ m (fcr)e"^ 
C 2 Z m+1 (kr)e^ m+1 ^ 



(21) 



where m is an integer, k > and Z m (x) can be any of 

the Bessel functions J m (x), Y m (x), and Hm (x). With 
the help of the well-known recursion relations of Bessel 
functions 



Z 'm( X ) ± — Z m{ x ) 
X 



one can show that the constants C\ and Ci satisfy 

/ k 2 2kk so \ ( d \ _ 2m*E 
\ 2kk so k 2 ) \ G 2 J ~ h 2 

Hence for a given k the two eigenenergies E± are 




E±(k) = 



2m* 



(k it fc so ) fc s , 



(22) 



(23) 



(24) 



Since the eigenvalues of the Schrodinger equation are in- 
dependent of the chosen coordinate systems the above 
eigenenergies should be the same as those given in 
Eq. (J2J, which is indeed the case when k = |k|. The 
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corresponding two non-trivial solutions for C^ 1 and 
are given by 

Cf/Cf = ±1, for E > 0, 



Cf/C* = -1, for - A so < E < 0. 



(25) 



For a given _E the two positive solutions of Eq. (|24[l for 
k are fc± given by Eq. (J2J. 

We are now in a position to construct different eigen- 
states using the Bessel and Hankel functions. The eigen- 
spinors regular at the origin are 



, r I Xm) = S 



±J m (k ± r) 
J m+1 (k ± r)e^ 

-J m (k±r) 
J m+1 (k±ryv 



e™ v , E > 0, 



e™ v , E < 0. 



(26) 



To derive the free-space Green's function in polar coor- 
dinates we shall also use solutions which are singular at 



the origin: 
(r|^) = 

(r\h+)= 
(r I h~) = 



±H£\k±r) 
H^ +1 (k ± r)e^ 



-H%\k+r) 



-Hil\k_r) 
H^ ] +1 {k^r)e^ 



E>0, (27a) 



', E < 0, (27b) 



E < 0. (27c) 



D. Free-space Green's function in polar coordinates 

Using and l|20b|) the free-space retarded Green's 
function in the two energy ranges becomes 



Goo(£;,r,r') = - 



Afi 2 k 



k-Hl + k+H\ 

liP (-§-r+ L r-k)t H -- H \) 

k-Hl + k+H 2 



k-Hi +k+H\ 



k-Hl +k+H\ 



E>0, 



-A so <E<0, 



(28) 



where we used the notations H± 2 = Hq (k±\r — r'|). 
In the off-diagonal elements the differentiations with re- 
spect to r and tp can be carried out by introducing a new 
variable p = r r'. Then, for E > our simple algebraic 
method yields the same result that was derived by Walls 
et al4i using a different approach. However, they do not 
present any explicit form for E < 0. 

In our previous paper— we used another form for the 
free-space Green's function (although as a lack of space 
it was not published there) in order to determine exactly 
the Green's function for circular Rashba billiards. In this 
approach differentiations with respect to r and ip in the 
off-diagonal elements were performed using the addition 
theorem of the Bessel functions^ 



H^(\r-r'\)= H^ 2 \r)J m (r>y m ^'\ r > r' 

m— — oo 

(29) 

and the recursion relations l|22|) . Then, for the free-space 
Green's function in polar coordinates for r > r' we obtain 
a rather compact form in terms of the spinors defined in 



Po| |77|): 

OO 

Goo (r, r')=cJ2 [k+ (r\h+) ( X + |r') + fc_ (r\h m ) ( Xm W) , 

m— — oo 

(30) 

where c = —im* /{Ah k). We shall use this form in sub- 
section 

III. CIRCULAR RASHBA BILLIARDS 

We now consider a circular Rashba billiard of radius 
R. The eigenstates of the system can be written as a 
linear combination of the regular eigenspinors given by 
Eq. i|26|) and the linear combination coefficients are cho- 
sen such that the eigenstates satisfy the Dirichlet bound- 
ary conditions. The straightforward calculation yields 
the following secular equation: 

J m {k + R)J m+ i{k-R)+sgn(E)J m (k-R)J m+1 (k+R) = 0, 

where m is an integer. For each quantum number m the 
solutions of this equation for E give the energy levels of 
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the circular Rashba billiards. The same secular equation 
was derived in Refs.i^SLS This equation is invariant 
under the change to — » —to — 1 (Kramers degeneracy). 
Formal solutions of the secular equation having zero wave 
vectors fc+ or fc_ are excluded since the corresponding 
wave functions vanish everywhere inside the billiard. One 
such solution is at E = — A so . 

Following the ideas of the systematic method of Berry 
and Howls' we have calculated the first few leading 
terms of the smooth counting function N{E). To do this 
we need the exact Green's function for circular Rashba 
billiards which is calculated in the following subsection. 



A. Green's function for circular Rashba billiards 

Boundary conditions for billiards requires that the 
Green's function should vanish at the boundary (i.e., if 
either r or r' is on the perimeter). The free-space Green's 
function ijHUI) for a given energy usually does not vanish at 
the boundary of the billiard. To fulfill the billiard bound- 
ary conditions we look for the exact Green operator, as 
usual, in the form of G = + Gh, where the homoge- 
neous Green's function satisfies [E — H ) Gh = 0. The 



boundary conditions for G are 



G(r, 



= G oc (r,r')+G H (r,r') = 0, for |r| = R, 

(32) 

where r' is inside the billiard. Since the homogeneous 
Green's function Gh satisfies the same Schrodinger equa- 
tion as the regular solutions given by Eq. (|26|l one can 
construct Gh from these eigenstates as 



Gh 



An 



Xm)(Xm I + B m \ Xm) (Xm I 



+Cm | Xm) (Xm I + D m I Xm) (Xm I 



(33) 



where the constants A m , B mi C m , and D m , in principle, 
can be determined from the boundary conditions (|32|1 . 
For arbitrary shapes of billiards it results in an infinite 
set of linear equations for the constants. Fortunately, 
in case of circular billiards the constants can be deter- 
mined analytically. Indeed, substituting Eqs. I|3UI) and 
(|33[l into Eq. I|32[) . and identifying the coefficients of the 
eigenspinors (x m | r') one finds 



A n 



Cn 



Xr, 



:k + {r\h m ), (34a) 



Xr, 



-C ki- 



ll; 



(34b) 



where the eigenspinors are evaluated at |r| = R. These 
equations, in fact, are four independent linear inhomo- 
geneous equations for A m , B m , C m , and D m since each 
eigenspinor is a two component vector. The solutions 
can be easily obtained from the appropriate determinants 
formed from the coefficients of Eq. I|34|) , and are given by 



4 — — - 



cfc+ 

± m 



D, 



c k- 



H { rk\k + R) -J m (k-R) 
H^ +1 {k+R) J m+ i(k-R) 

-H%\k+R) -J m (k-R) 
H ( V +1 (k + R) J m+ i{k-R) 



J m (k+R) -H&\k-R) 

J m+1 (k+R) H^ik-R) 

-J m (k+R) -H£\k-R) 

J m+1 (k + R) H { V(k-R) 



B n 



2ic 

TTRF m ' 



where 



E > 0, 

,E<0, 

(35a) 

,E>0, 

,E<0, 

(35b) 

(35c) 



F m = J m (k+R)J m+1 (k-R)+sgn(E)J m (k-R)J m+1 (k+R). 

(35d) 

In Eq. H35cfl w e used the Wronskian relations for the 
Bessel functions^ 

Finally, the analytical form of the exact retarded 
Green's function of circular Rashba billiards in polar 
coordinates is a sum of the free-space Green's function 
Goo(r, r') given by Eq. pgjl or Eq. lj3U|) . and the homoge- 
neous part G^(r, r') = (r | Gh | r'), where the operator 
Gh is given by Eq. $£ft together with Eq. JSgJl. 

The eigenenergies of any billiards can be obtained from 
the poles of the retarded Green's function G. For circular 
Rashba billiards the poles of G are the poles of Gh , i.e., 
the zeros of F m . As can be seen it yields the same secular 
equation (|31|) derived independently, and thus it provides 
one check point for the Green's function Gh- 



B. The smooth part of the density of states 

To calculate the DOS and the counting function for 
circular Rashba billiards we adopt the ideas of the sys- 
tematic method of Berry and Howls.— The exact Green 
operator of the system is G — G^ + Gh ■ The first term 
of the density of states © is the contribution from Goo 
in the trace of G. The result is given in l|14|) . while the 
leading term in the counting function N(E) is given by 

03- 

The correction terms of the DOS can be obtained from 
the trace of Gh given by Eq. (|33|l . This involves the limit 
r — > r', the trace of the 2 by 2 matrix Gh (r, r) (trace over 
spinor indices) and the integration over the area of the 
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billiard. After a straightforward calculation we found 



^ pit 
Tr Gh = 27r / rdr 

m _ r*r>, « 



m— — oo 



+J^(fc_r)(D m _i + An) 



+2J m (fc+r)J m (fc_r) 



\B m -i — B m ), E > 
(B m -i + B m ), £<0 



(36) 



In the series with terms A TO _i the summation index m 
has been shifted by one to have the same radial integral 
as that in the series for A mi and the same trick was done 
for series containing B m -i, C m -\ and D m -i. 

The radial integrals in H36|) can be performed analyt- 
ically. 4 — To calculate the density of states one needs to 
evaluate Tr Gh at complex energies E + ir\. To this end 
we follow the approach originally applied by Stewart- 
son and Waechter^i and later for example Berry and 
Howls j2£ and the Bessel functions of the first kind J m {z) 

( 1 2) 

and H^\z) are converted to the modified Bessel func- 



tions I m (z) and K m {z) by extending the energy E to the 
complex plane. This is the so-called heat-kernel method. 
However, in our case one has to be careful for negative 
energies. It turns out that the parameters x, x+ and X- 
dcpcnding on energy E (here E is real) and defined as 



ix = Rk(E + ir]), 
ix + = sgn(E)Rk + (E + ii]), 



IX- 



Rk-{E + irj), 



(37a) 
(37b) 
(37c) 



are useful to convert the Bessel functions of the first kind 
to the modified Bessel functions using the identities 

J m (iz) = i m I m {z), (38a) 
if« (iz) = - (-i) n+1 K m (z) , -Tr < arg z < ~ (38b) 

7T Z 

H%\-iz) = -i n+1 K m {z), 4<argz<.. (38c) 

7T Z 

After a tedious algebra the result of these transformations 
in Eq. (|36|l can be written as 



TrG H {E + i V ) = 1 ^- 
n x 



f m (E + ir)), where 



m— — oo 



fm(E + irf) 



m 



x + I m (x+)K m (x + ) 



m 



x-I m (x-)K m (x-) 



-Pjn (x~\- ) X— ) 



m 



I m {x+) 



m 




I L( x +) _ x -CQe-) 

I m (x+) ~ I m (x-) 



,(39a) 



and 

P m (x+, X-) 



1 



I 'm(x+) + I' m {x-) _ m__ m_ + I' m (x + ) + I' m {x-) 
I m {x+) I m (x-) x + x- 



r 



rn 

x+ 



rn 

X- 



(39b) 



Up to now the trace of Gh for circular Rashba billiards 
is exact. Note that the above mentioned transforma- 
tions result in the same form of Tr Gh (E + irf) both 
for negative and positive energy E. Moreover, as a self- 
consistent check, one can show that f m (E) becomes the 
same as that in Refs.^2S when the spin-orbit coupling 
is zero, i.e., for x+ — > X-. Indeed, in this limit, using 
the L'Hospital's rule and the Bessel differential equation 
for I mi it can be shown that the factor multiplied by 
P m (x+, X-) in Eq. I)39|) is exactly zero, and the remain- 
ing terms can be rewritten in the same form as that in 
Refs.ii^ 

The next step is to replace the modified Bessel func- 



tions in Eq. (|39|l by their uniform approximation^^ 
Keeping only the leading terms we obtain 



m*R 2 



TtG 



H 



E 



m— — oo 



l 



X+ + X- 



1 - 



l^J(m 2 + x 2 + ){m 2 + x 2 _) 



(40) 



Note that the second term bracketed in square brackets 
is zero when x+ — > X—, and one finds the perimeter term 
of the DOS for billiards with zero spin-orbit coupling 
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from the remaining terms. Taking into account the sub- 
sequent terms in the uniform approximation provides a 
systematic way to derive higher order terms for the trace 
of Gh as in Refl32l for normal circular billiards. How- 
ever, with non-zero spin-orbit coupling the calculations 
become more cumbersome. 

The summation over m in l|4(J|) can be rewritten using 
the Poisson summation formula 28-44 



E 



/«*= E 



dmf n 



(41) 



Then, the Weyl series, i.e., the smooth part of the DOS, 
following Berry and Howls^ can be obtained by keeping 
only the [i = 0-term in (|41|l . Carrying out the limiting 
process, rj — > in the trace of Gh given by Eq. H4()(l . and 
using the integral 



dz ■ 



2(a + b) 



2z 2 -a 2 -b 2 
^/(z 2 -b 2 ){a 2 ^ 2 ) 
' a — 6 N 



E 



a + b 



- K 



a + b 



(42) 



valid for < b < a (E and K are the complete elliptic 



integrals with the same definitions as in Ref|43J) for inte- 
gration over m, we obtain the contribution to the smooth 
DOS coming from Ti Gh- A tedious calculation yields 



Qh{s) 



1 



e+e s . 



2^/e + e so 



E 



e+e 



, £ > 0, 

(43) 

, e < 0, 



where the dimensionless energies e = 2m*ER 2 /h 2 and 
s so = 2m*A so R 2 /h 2 — k 2 Q R 2 have been introduced. The 
first term is the contribution from the first and second 
terms of Eq. (|40|l . and it coincides with the perimeter 
term derived in Eq. 118|) for arbitrary shapes of Rashba 
billiards. The Dirac delta term and the terms containing 
the complete elliptic integrals in (|43|l come from the term 
involving bracets in (|40|) . 

Finally, including the contribution from Tr Goo , the in- 
tegration of the DOS over E leads to the smooth counting 
function N(e) for Rashba billiards: 



£+2e ao 
2 



N(e) 



EL/^±— ) , for -e so <e<0 



, for £ > 0, 



(44) 



The first two terms (for both positive and negative ener- 
gies) are the contribution from Goo. They are the area 
and perimeter terms in the Weyl series and agree with the 
results given by (fL5j) and (|T^1) . respectively for arbitrary 
shapes of Rashba billiards. The terms containing the 
complete elliptic integrals are corrections to the perime- 
ter term in case of circular billiards. We note that in 
a completely different context, namely for annular ray- 
splitting billiards, a similar Weyl formula has been cal- 
culated^ involving also elliptic integrals. 

We have compared the smooth counting function N(e) 
given by Eq. (|44|l with the exact counting function N(s) 
calculated from the energy levels obtained from the sec- 
ular equation l|31|l for different m. The relevant parame- 
ter characterizing a circular Rashba billiard of size R is 
k so R- Typical values for the spin-orbit-induced spin pre- 
cession length L so = n/k so are of the order of a few hun- 
dred nanometers i Taking R — 10/im for a typical size 
of quantum dots, the relevant parameter k so R in Rashba 
billiards can be as large as 70. Figure^, shows the exact 
and the smooth counting functions as functions of the 
dimensionless energy e. There are 6388 energy levels in 



the plotted energy range. To see the difference between 
the two functions, in the inset we plotted them close to 
the bottom of the energy spectrum. Figure^) shows the 
difference AN — N(e) — N(e) as a function of e. The 
difference fluctuates around zero, which means we did 
not miss levels (the mean value of AN is a sensitive test 
for missing levels, see e.g., 

ReflU). 

Without correction 
terms in Eq. I|44|l with elliptic integrals, A7V would in- 
crease monotonically on average as shown in Fig. [21 and 
would predict a difference ~ 27 in the energy range plot- 
ted. 



C. The counting function for negative energies 

In Fig- El the exact counting function is shown for neg- 
ative energies near the bottom of the spectrum — £ so . As 
can be seen the exact N(e) shows an additional rounded 
step structure at certain energies £* . This feature shows 
up only for negative energies, although for larger ener- 
gies this is less pronounced. The step structure results 
in large deviations A7V at energies £* and concomitant 
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FIG. 1: (Color online) In panel a) the exact counting function 
N(e) (dashed line) and N(e) (solid line) are shown for y 7 ?^ = 
k ao R = 70. The inset shows the enlarged portion of the main 
figure close to the bottom of the spectrum. In panel b) the 
difference A7V = N(e) — N(e) is plotted. In both panels 
dimensionless energies e = 2m*ER 2 /h 2 are used. 




3000 



FIG. 2: The difference AN between the exact counting 
function and N(e) without the terms containing the ellip- 
tic integrals in Eq. H44|l is plotted. The energy is scaled as 
e = 2m*ER 2 /Ti 2 and k EO R = 70. 



large peaks in the DOS. 

To see the reason for this behavior, it is useful to plot 
the energy levels as functions of m, as shown in Fig. 0] 
The curves in the figure start almost horizontally at £* , 
n = 1,2,... resulting in large peaks in the DOS at the 
same energies. Using Debye's asymptotic expression for 
Bessel functions with large argument^ we were able to 
derive the energy dispersion in leading order: 



(ff 



(45) 



valid only for negative energies. Figure 01 also shows 
the comparison of the exact energy levels and their ap- 
proximated m and n dependence given by Eq. 145fl . For 
small 77i, 7i the above expression agrees excellently with 
the numerics (e.g., £o,i is accurate up to 7 digits for 
e so = 70). The smallest energy level in the spectrum of 
the circular Rashba billiard is E m \ n — h 2 /(2m*R 2 ) e ,i — 
?7 2 /(2m*)7r 2 /(4i? 2 )- A so . 

We now derive an approximated expression for the 



FIG. 3: (Color online) The exact counting function N(e) 
(solid line) and the asymptotic counting function A^aaymp(e) 
given by Eq. (gHJ (dashed line). The energy is scaled as 
e = 2m*ER 2 /h 2 and k BO R = 70. 



-4700 



-4750 



£ -4800 



-4850 - 
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FIG. 4: (Color online) The m dependence of the exact en- 
ergy levels (in units of Tl 2 /2m* R 2 ) of circular Rashba billiards 
(symbols) for a given n ranging from n = 1 to n = 7. The 
solid lines are the curves obtained from the approximation of 
the exact energy levels given by Eq. I|45|l as functions of m 
with the corresponding n. Here k so R — 70. 



counting function using 



iVa. 



>(e) 



E E 

m— n—1 



0(e 



(46) 



where £ m , n are given by Eq. 1)450. the factor 2 takes 
into account the Kramers degeneracy in m, and m max — 
[•y^so ] and n max = [(2y / £^/7r] are the largest m and n for 
which e m . n is still negative. Here [•] stands for the inte- 
ger part. Applying the Poisson summation formula™^ 
in the sum over m in Eq. (|46|l and keeping only the non- 
oscillating term we find 



,(e) = 2 E / x Q{e - e m . n )dm 

n=l " 2 
n max 

= 2 m*{s,n), (47) 
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where m*(e,n) is the solution of e m>n — e for m at a 
given e and n. Thus, from Eq. i|47p. after some sim- 
ple algebra, we obtain the final form of the asymptotic 
counting function in Debye's approximation: 



c± satisfy 



^asymp(^) — 2 £ so 



where et 



E 

n=0 

m 2 



■Q(e 



for e < 0, 



(48) 

The result is plotted 



together with the exact counting function in Fig. The 
agreement is clearly visible near the bottom of the spec- 
trum. However, it is an open question what semi-classical 
picture can be associated to the content of Eq. I|45[) . A 
possible treatment in this direction may be the semi- 
classical approach of ReflZoL 

The density of states is the derivative of the counting 
function N(E) with respect to E, therefore for circu- 
lar Rashba billiards in the DOS square root types sin- 
gularities (van Hove type) appear at energies E^ ng = 
S 2 /(2m*i? 2 )e*. Attaching leads to a circular Rashba 
billiard, the transport properties of this open system are 
determined by the tunneling conductance which is pro- 
portional to the DOS. Thus, the measured conductance 
should be changed abruptly in the negative energy spec- 
trum of circular Rashba billiards at energies E^ ns . 



D. The spin structures of the eigenstates 

It is straightforward to obtain corresponding spinor 
eigenstates and calculate their expectation value for the 
z component of spin. Similar to the case of Rashba-split 
eigenstates in rings^ but in contrast to that of quantum 
wiresji2*ii it turns out to be finite. 

The eigenstates of the Rashba billiards satisfying the 
Dirichlet boundary conditions can be expressed with the 
linear combination of the regular eigenspinors | Xm) given 
by HJ: 



J m (k-r) \ I imip 
' 1 J m+1 (k-r)e^ \ e ' 



(49) 



J m (fc_i?) J m+1 {k_R) 
J m {k+R) J m+ i(k + R)' 



(50) 



and k± satisfy the secular equation l|31[l with energy lev- 
els e m ,n- Eigenstates given by Eq. (|4*9"|l are valid for 
£ m ,n > 0. In the opposite case, one should use the 
eigenspinor given in l|2t)|l for E < 0. Regarding the spin 
structures, it turns out that both cases (the positive and 
negative energy levels) can be treated at the same level if 
the definitions for k± in J3J arc modified as k± = k^fk so . 
Therefore, hereafter we use these new definitions for k±. 

The spin structure of the eigenstates in Rashba bil- 
liards can be obtained by calculating the expectation val- 
ues for spin components: 



(*<), 



R r 2n 



rdrdtp*+ in {r,<p)ai* m>n (r,<p), (51) 



where i = x,y,z, and + denotes the transpose and the 
complex conjugation of a spinor state. The integrand in 
this equation is the spin density of Cj. The eigenstates 
(|49|l can be written in the form of 



imtp 



(52) 



where Af is the normalization constant, the coefficients 



and then it is easy to show that the spin density depends 
only on r as 

*+ n (r)a,* miB (r) = |*«„(r)| 2 - |*£> B (r)|*(53a) 

*+ B (r)«7 P * m , n (r) = 2*£> n (r)*« n (r), (53b) 

where ay = cos tpa x + sin tpa y is the in-plane radial com- 
ponent of the spin. One can also show that the angu- 
lar component \P+ n a v ^ m ,n of the in-plane spin den- 
sity is exactly zero, where a v — — simp a x + cos ipa y . 
Therefore, the in-plane spin density at point r in the 
billiard is along the radial direction This implies 
that the expectation values for in-plane spin is zero, i.e., 

(^*x)m,n — (^y)m.n 0. 

Performing the integration (that can be carried out 
analytically) in Eq. i|51[l for (u z } m ,n we find 



Jm(k-R) 

J m+ x{k-R) 



J m + 1 (k-R) 
J m (k-R) 



(54) 



This is an exact analytic result for the expectation val- billiards. 

ues of the z component of the spin for circular Rashba Figure shows the expectation values of (a z ) m cal- 
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culated numerically from (|54|) for different angular quan- 
tum number m and eigenvalues e mi „ with a given Rashba 
coupling strength a. One can see from the figure that 




FIG. 5: The expectation values of (a z ) n as functions of the 
angular quantum number m and eigenvalues e m ,n f° r k so R — 
70. 

(o~ z ) m n has a peak at m — and for eigenvalue e mj „ close 
to zero. We have studied how this peak value changes for 
different Rashba coupling strength a. For each k so R and 
m the maximum of (o~z) m n over the eigenvalues e m ,n is 
plotted in Fig. It is clear from the figure that the ex- 
pectation values of (o~z) m is robust for different Rashba 
coupling strength a. 



1 
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0.6 
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man effect. Using the first order perturbation valid for 
weak field limit, i.e., when the cyclotron radius is much 
larger than the size of the Rashba billiard, the values of 
the Zeeman splitting is proportional to the expectation 
values of \Cz) • Thus, we believe that the significant 
magnitude of the spin z component found from our nu- 
merical results can be detectable experimentally. 



IV. CONCLUSIONS 

Before concluding in this section we briefly summarize 
our results not discussed in this paper on the statistics of 
energy levels, and highlight some open theoretical prob- 
lems in connection with Rashba billiards. 

The Schrodinger equation (including boundary condi- 
tions) for circular Rashba billiards is separable in po- 
lar coordinates, thus integrable. Hence, the statistics 
of energy levels should be Poissonian (see e.g. Ref®. 
Indeed, we have found that the nearest-neighbor level- 
spacing distribution P(s) is Poissonian (not shown here). 
For other shapes of Rashba billiards, spin-orbit coupling 
may destroy integrability, in which case Random Matrix 
Theory (RMT) predicts that the level statistics should 
be governed by the symplectic ensemblei244I Note, how- 
ever, that some intermediate distribution (not described 
by RMT) was found^ for a rectangularly shaped billiard 
in the limit of small k so , reflecting the fact that a rectan- 
gular billiard is integrable in the absence of SO coupling 
but non-integrable when SO is finite. 

We now list a few interesting open theoretical prob- 
lems. The Weyl formula is essential to develop a periodic 
orbit theory for Rashba billiards. (For normal billiards, 
see Brack and Bhaduri's book in ReflHH and a the- 
ory in case of harmonically confined Rashba systems is 
given in Ref 40.) The Green's function method presented 
in this work would be a suitable starting point to calcu- 
late observables such as the magnetization^ or persistent 
currents^ in Rashba billiards. 

In conclusion, we have presented a study of electron 
billiards with spin-dependent dynamics due to Rashba 
spin splitting. Semi-classical results for the spectrum 
agree well with exact quantum calculations. We find in- 
teresting properties of negative-energy states, including 
a finite spin projection in the out-of-plane direction. 



FIG. 6: The maximum of the (a z ) as functions of m and 

\ / m , n 

For Rashba billiards, in weak magnetic field the energy 
levels of the Kramers doublets will be splitted by the Zee- 
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